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CO ' Abstract. We present new photometric data for seven high-amplitude d Scuti stars. The observations were acquired between 

, 1996 and 2002, mostly in the Johnson photometric system. For one star (GW UMa), our observations are the first since the 

discovery of its pulsational nature from the Hipparcos data.The primary goal of this project was to update our knowledge on the 
period variations of the target stars. For this, we have collected all available photometric observations from the literature and 
^ . constructed decades-long 0-C diagrams of the stars. This traditional method is useful because of the single-periodic nature of 

the light variations. Text-book examples of slow period evolution (XX Cyg, DY Her, DY Peg) and cyclic period changes due 
^ , to light-time effect (LITE) in a binary system (SZ Lyn) are updated with the new observations. For YZ Boo, we find a period 

5_( ' decrease instead of increase. The previously suggested LITE-solution of BE Lyn (Kiss & Szatmary 1995) is not supported with 

^ ' the new O-C diagram. Instead of that, we suspect the presence of transient light curve shape variations mimicking small period 

, changes. 
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1 . Introduction In Paper I (Kiss et al. 2002b) we discussed the double-mode 

pulsation of V567 Ophiuchi. The aims of this project have been 
High-amplitude 5 Scuti stars (hereafter HADSs) are either Pop. .,• j • r^- ^ i /onmuN j» a e 

° ^ ^ ^ F outlined in Kiss et al. (2002b) and to avoid repetition, we refer 

I stars close to the main sequence or evolved Pop. II stars (these , i j u i . j- j • »u . 

^ ^ ^ to theoretical and observational aspects discussed in that pa- 

are also known as SX Phe variables) with very characteristic rru ■ ■ f.u- • . uv u J . 

' •' per. The main aim of this paper is to publish a new dataset for 

light variations caused by radial pulsations (Rodriguez et al. • ui . n j u . . • j . 

° ^ ^ seven variable stars. Besides the new photometnc data we give 

1996). Typical periods range from 0t05 to 0?15 associated with a ^ ■ e^u ■ ■ a u ■ u u • -ru 

^ f & an updated description of their period changing behaviour. The 

relatively large amplitudes (Av > OT'SO is a widely adopted ^ ■ t a ^ ■ ^ ^ u • i 

J ^ f \v- J f analysis of multicolour data in terms of physical parameters 

convention). Owing to their short periods and high amplitudes, .,, , uv u j • ^ i c^i. ■ r-r^ t . i 

^ ° tor-' will be published in the final paper of the series (Derekas et al., 

these stars are very good targets for small and moderate-sized . . 

° ° in prep.), 

telescopes, so that interesting astrophysical phenomena can be t,. ^ . ^ • * n • 4. i 4.- 

^ to r- ^ r 'pjjg structure of the paper is the following. Target selection 

easily studied even with modest instrumentation. Apart from , , , u j • c » o t-u • ^ ci, 

■' and observations are described in Sect. 2. The main part of the 
attempts to detect evolutionary effects, the most interesting u-u.u 1. .jc 

^ -' ° paper is Sect. 3, in which the results are presented. Separate 

case studies are related to the suspected binary HADSs, in u j . j . • j- j 1 * u 1 

-' subsections are devoted to individual stars, where every rel- 
which the physical parameters of the pulsating component can ^ j . 1 / 1. ^- 1 j ^ 1 v j . j 

^ ^ ±- o ±- evant detail (observational data, sample light curves, updated 

be constrained from the binary nature. However, the time-span r\ a- \- ■ Au f • .j c. 

-' ^ O-C diagrams) is given. A brief summary is presented in Sect, 

of the available observations is often too short, therefore, sig- ^ 

nificant observational efforts have to be done to reach unam- 



biguous conclusions. 



2. Target selection and instrumentation 



Send offprint requests to: L. L. Kiss, 
e-mail: laszloOphysics. usyd.edu.au The observations were started with regular monitoring of 

* On leave from University of Szeged, Hungary BE Lyncis, for which we suspected binarity (Kiss & Szatmary 
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Table 1. The list of programme stars. The magnitude ranges 
and periods are taken from the GCVS. 



Star 


Pop. 


^max 




P(d) 


Type of obs. 


DYHer 


I 


W15 


10?66 


0.14863 


V 


YZBoo 


I 


IGP30 


icrso 


0.10409 


V 


XXCyg 


II 


ir28 


12? 13 


0.13486 


V{RI)c 


SZLyn 


I 


9.''08 


9^12 


0.12053 


V 


BELyn 


I 


8°60 


9?00 


0.09587 


BV/uvby 


DYPeg 


n 


9?95 


10P62 


0.07293 


BV/uvby 


GWUMa 


n? 


9?48 


9^91 


0.20319° 


BVIc 



" ESA (1997) 



1995). The first extension was made toward XX Cygni (Kiss 
& Derekas 2000) and since then, we have tried to include all 
HADSs in the northern sky brighter than V = 1 1™3 at max- 
imum (this Umit was determined by the typical accuracy and 
limiting magnitudes of our main instruments). Unfortunately, 
this goal was unreachable because of unfavourable weather 
conditions in several observing seasons. That is why we 
excluded the following bright northern HADSs: GP And, 
AD CMi, DH Peg, V1719 Cyg. Also, we concentrated on the 
monoperiodic variables and some of these excluded stars are 
well-known double-mode pulsators (see Rodriguez et al. 1996). 
Thus, our sample consisted of 10 variables, however, two stars 
(V 1 1 62 Ori and DE Lac) had so meagre coverage that we had to 
remove them from the final target Ust. V567 Oph turned to be a 
double-mode pulsator and was discussed in Kiss et al. (2002b). 
The remaining seven HADSs are listed in Table 1, where their 
main observational properties are also summarized. 

During the seven years of observations, we have uti- 
lized various instrumentations at three observatories (Szeged 
Observatory, Piszkesteto Station of the Konkoly Observatory, 
Sierra Nevada Observatory). In the following we briefly intro- 
duce the telescopes and detectors used in this project. As the 
primary aim was to obtain good light curve coverage enabling 
accurate determination of the epochs of maximum, the CCD 
measurements were uninterrupted single-filtered V-band obser- 
vations, with a few exceptions. 

- Szeged Observatory, 0.28m Schmidt-Cassegrain(Sz28) 
This telescope is located in the very centre of the city of 
Szeged, thus suffering from strong Ught pollution. For our 
observations on two nights in 2002, it was equipped with 
an SBIG ST-71 CCD camera (765x510 9yum pixels giving 
a field of view (FOV) of li;5x7.'5). This instrument was 
the least used in our project. 

- Szeged Observatory, 0.4ni Cassegrain (Sz40) 

The majority of the observations were acquired with the 
0.4m Cassegrain-telescope of the Szeged Observatory. 
Between 1996 and 2000, photoelectric photometry was car- 
ried out through Johnson BV filters. The detector was an 
SSP-5A photoelectric photometer and we obtained differ- 
ential photometric data using selected comparison stars lo- 
cated near to target stars. The photometer was replaced by 
an SBIG ST-9E CCD camera in 2001 (5 12x5 12 20//m pix- 
els, FOV=6' X 6') and most single-filtered V-band observa- 
tions were made with this instrument. 



- Piszkesteto, 0.6m Schmidt (P60) 

Fewer observations were done with the 60/90/ 180cm 
Schmidt-telescope mounted at the Piszkesteto Station of 
the Konkoly Observatory. The detector was a Photometries 
AT200 CCD camera (1536x1024 9//m pixels, FOV= 28' x 
19'). The observations were done either in B or V filters. 

- Sierra Nevada Observatory, 0.9m Ritchey-Chretien 
(SNOW) 

Occasionally, we acquired simultaneous uvby Stromgren 
photometric observations using the 0.9m telescope of the 
Sierra Nevada Observatory (Spain) equipped with a six- 
channel Stromgren-Crawford spectrograph photometer. 

The data were reduced in a standard fashion. For the photo- 
electric observations, we made use of different computer pro- 
grams written by LLK, which were based on methods de- 
scribed in Henden & Kaitchuk (1982). The CCD observations 
were reduced with standard tasks in IRAF', including bias re- 
moval and flat-field correction utilizing sky-flat images taken 
during the evening or morning twillight. Differential magni- 
tudes were calculated with aperture photometry using two com- 
parison stars of similar brightnesses. We have not tried PSF 
photometry because the small FOVs contained too few stars 
for a reliable PSF-determination (see Vinko et al. 2003 for a re- 
lated discussion using the same instruments). The typical pho- 
tometric accuracy was +0^01 -0^02, depending on the weather 
conditions and brightness difference of the variable and com- 
parison stars. In fields far from the galactic plane - e.g. those 
of BE Lyn and SZ Lyn - we had to choose somewhat fainter 
comparisons, typically 2-3 mags dimmer than the target, which 
consequently increased the noise level of the data. 

Thanks to the dense sampling of the light curves, new times 
of maximum light were easy to determine from the individual 
cycles. This was done by fitting low-order (3-5) polynomials 
to the light curves around maxima. We repeated the fitting pro- 
cedure by changing the parameters of the fits to yield some in- 
sights into the uncertainty range of the determined epochs. We 
estimate the typical accuracy as +0?0003 (i.e. +26 sec), which 
is comparable with the usual exposure lengths. 

3. Results 

The period variations were studied with the classical O-C (ob- 
served minus calculated) method. The use of this simple ap- 
proach is justified by the remarkable stability of the HADS 
light curves, both in amplitude (Rodriguez 1999) and in period 
(Rodriguez et al. 1995, Breger & Pamyatnykh 1998, hereafter 
BP98). The blind use of the O-C technique has been strongly 
criticised by Lombard and Koen in their series of papers (e.g. 
Lombard & Koen 1993, Koen & Lombard 1995, Koen 1996, 
Lombard 1998) underlining the importance of the evaluation of 
the model residuals. For this, they developed various statistical 
methods by analysing autocorrelation of the residuals. This is 

' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 2. Journal of observations and new times of maximum 
for DY Her 



Date 


filter 


Inst. 


Points 




2001-05-09 


V 


Sz40 


509 




2001-05-10 


V 


Sz40 


663 


2452040.3804 










2452040.5317 


2001-05-19 


V 


Sz40 


409 


2452049.4471 


2001-06-25 


V 


Sz40 


506 


2452086.4579 



-1 8' ^ ' 1 ' ^ ' ^ ' 1 

2049.4 2049.45 2049.5 2049.55 2049.6 

HJD (2450000+) 

Fig. 1. Sample light curve for DY Her. 



a crucial point when studying long-period stars (e.g. Mira vari- 
ables) with relatively large stochastic scatter of the period. In 
case of the HADSs one does not expect such intrinsic period 
scatter (as suggested by the observations) and evolutionary pe- 
riod changes might be constrained with the O-C technique. 

Nevertheless, in every case we have checked the model 
residuals via their frequency analysis. This can be used to 
search for secondary periodicity or slow changes of the period 
hidden by the observational scatter. The basic idea is the fact 
that a small modulation of the O-C diagram appears as an ad- 
ditional periodicity of the O-C residuals. A recent example for 
this has been given by Pocs et al. (2002), who presented the 
Fourier spectrum of the O-C diagram of the multiply periodic 
HADS RV Ari, clearly revealing fm - fx - fo in their data. 
On the other hand, if there is a slight period change in the data 
and a wrong model is subtracted, then the residuals' spectrum 
will show a peak very close to / = ^. This closely separated 
peak in the Fourier spectrum, however, can be explained by an- 
other way, too: it appears simply because of the periodicity of 
the spectrum itself. The Fourier transform of an O-C diagram 
is periodic with the frequency of the ephemeris because of the 
periodicity of data sampling. Therefore, if a long-term trend is 
present in the residuals, that causes a low-frequency component 
at fiovi, which appears again at /low + /• Whatever explanation 
we adopt, Fourier analysis of the residuals can help emphasize 
the hardly visible trends of the data. 



0.04 




50000 

cycle number 



Fig. 2. The O-C diagram of DY Her. The solid line shows the 
parabolic fit. 



were used (comp^GSC 0968-1532, check^GSC 0968-1002, 
both stars are about 1 1™0). The full log is given in Table |2 
a sample light curve is shown in Fig.^ There is no new time 
of maximum for the first night, because we had to interrupt our 
observations for half an hour and the maximum just appeared 
during that break. 

We have collected all times of maximum from the literature 
(the updated list of references is given in Pocs & Szeidl 2000) 
to arrive to the final O-C diagram. It is plotted in Fig. |2 and 
was calculated with the following ephemeris: 



3.1. DYHerculis 

The most recent period-change study for DY Her was made by 
Yang et al. (1993), but even the last data of that paper were ob- 
tained in 1984. Previously, Szeidl & Mahdy (1981) collected 
all available data to study long-term period variations. Both 
studies revealed a slow period decrease of the star. The main 
characteristics are summarized by Pena et al. (1999). Recently, 
Pocs & Szeidl (2000) interpreted the long-term O-C diagram 
in terms of the light-time effect (LITE) caused by a hypotheti- 
cal low-mass component, although their conclusions were quite 
ambiguous. 

We observed the star on four nights in 2001 with Sz40. 
More than 2000 individual exposures of 20-40 seconds were 
obtained in V-band. Two comparison stars located within 6' 



i/yD^ax = 2433439.4871 + 0.1486309 x E 

(the epoch is from the GCVS, the period from ESA 1997). 
The parabolic fit has the form 

O-C ^ -0.0001(2) + 3.9(1) X 10"^£ - 8.3(1 1) x 10'^^ 

with an rms of 0?00115 (ierrors of the last digits are in 
parentheses). The linear fit yields higher rms (O?00157) and 
the Fourier spectrum of the residuals (Fig. |3} shows that the 
parabolic residuals are closer to an uncorrected noise. The 
peak in the spectrum of the linear residuals appears close 
to the pulsational frequency (with wich the O-C was calcu- 
lated): the difference is 2.30 x lO"'* d^'. Therefore, the hypo- 
thetical binary nature (Pocs & Szeidl 2000) of the star is not 
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freq. (c/d) 



Fig. 3. Top panel: the Fourier spectrum of the residuals for 
DY Her after subtracting a linear fit. Bottom panel: the same 
with the parabolic fit. 



confirmed and we are confident about the definite detection 
of the slow period decrease. The parabolic coefficient (usu- 
ally denoted as /3/2-0.5 P dP/dt) results in a period change 
- (-2.8 + 0.4) X 10"^ year"', which is slightly smaller and 
more accurate than the value used by BP98 (-4 x 10 year '). 
Nevertheless, the main conclusion drawn by several earlier au- 
thors is not changed, as this smaller decreasing rate is generally 
in good agreement with the evolutionary calculations. 



3.2. YZBootis 

YZ Boo has quite a long observational record and the most 
relevant information on the star was summarized by Peiia et 
al. (1999). Its period change was studied by Szeidl & Mahdy 
(1981), Jiang (1985) and Hamdy et al. (1986). The latter au- 
thors analysed the longest dataset available at that time and 
concluded that the O-C diagram can be fitted with a positive 
parabola, indicating a constant period increase. They arrived at 
a relative period changing rate of = H-3 x 10"*^ year"'. 

Our observations were carried out using Sz40 on six nights 
in 2001 (Table The exposure time was between 15 and 40 
s. Unfortunately, only one comparison star was available in the 
small FOV of the instrument (GSC 2569-1 184 with V *1 1"3). 
Slightly more than 2500 data points were obtained. The longest 
subset, containing 697 points, is shown in Fig.|3 



Table 3. Journal of observations and new times of maximum 
for YZ Boo 



Date 


filter 


Inst. 


Points 




2001-03-16 


V 


Sz40 


279 


2451985.6116 


2001-04-01 


V 


Sz40 


536 


2452001.5381 


2001-04-25 


V 


Sz40 


177 


2452025.3695 


2001-04-28 


V 


Sz40 


697 


2452028.3902 










2452028.4946 










2452028.5975 


2001-04-30 


V 


Sz40 


526 


2452030.3693 










2452030.4726 










2452030.5759 


2001-05-09 


V 


Sz40 


307 


2452039.3221 



-1.2 



> 
<J 



YZBoo 2001.04.28/29. 



A 



\ 



V 



1. ' 

V 



2028.3 



2028.4 2028.5 
HJD (2450000+) 



Fig. 4. Sample light curve for YZ Boo. 
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Fig. 5. The O-C diagram of YZ Boo. The solid line denotes 
the parabolic fit. Residuals are shown below. 



Ten new times of maximum were determined from the indi- 
vidual light curves (TableO, which extended the existing data 
by another ~15 years. Published times of maximum were col- 
lected from Hamdy et al. (1986) and references therein, Kim 
& Joner (1994), Agerer et al. (1999) and Agerer & Hubscher 
(2000). The O-C diagram was calculated with the following 
ephemeris: 
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Fig. 6. Top panel: Fourier spectrum of the residuals for YZ Boo 
after subtracting a linear fit. Bottom panel: the same with the 
parabolic fit. 



//7Dmax = 2438206.6961 + 0.1040920 x E 

and is plotted in Fig. |5] The solid hne corresponds to the 
parabolic fit: 



O - C = 0.0018(2) - 4.25(6) x IQ-'^E - 1.5(6) x lO"'^^^ 

with an rms of O?00169. The data can be fitted almost equally 
well with a simple linear term (rms=O?00173), but a frequency 
analysis of the residuals slightly favours continuous period 
change. It is shown in Fig.|6l There is a peak at / ^ 9.6116 
d"'(~ -p + 4.7 x lO""* d"') in the spectrum of the residuals of 
the Unear fit, while the spectrum of the parabolic residuals are 
more noise-like. This suggest a real period change over the six 
decades of observations, although the evidence is fairly weak. 
If we accept its significance, the given second-order coef- 



ficient corresponds to a relative period changing rate - 
(-1 + 0.4) X 10"^ year"', which has an absolute value of about 
one third of the previous estimate by Hamdy et al. (1986) and 
slightly larger than that of Peniche et al. (1985), but with the 
opposite sign. It means that contrary to Hamdy et al. (1986) 
and Peniche et al. (1985), we found period decrease instead 
of increase. The determined rate has great uncertainty and we 
recall a note by Peniche et al. (1985): "...it is impossible right 



Table 4. Journal of observations and new times of maximum 
forXXCyg 



Date 


filter 


Inst. 


Points 


Anax 


2002-07-20 


V 


Sz40 


60 


2452476.4983 


2002-07-20 


Rc 


Sz40 


61 


2452476.4982 


2002-07-20 


Ic 


Sz40 


59 


2452476.4976 



0.5 



1.5 



AI (mag) 
AR (mag 
AV (mag) 



2476.48 2476.52 2476.56 

HJD (245000+) 



2476.6 



Fig. 7. V{RI)c light curves of XX Cyg. 



now to decide if the period is constant or if it is varying, but this 
will be feasible only if the star is regularly observed during the 
next 40 years". Almost 20 years have passed since then and the 
uncertainty is still too large. That clearly shows the need for 
continuous follow-up in the future. The given slow rate is in 
good agreement with theoretical calculations of BP98, which 
predict small period changes of the same order of magnitude. 



3.3. XX Cygni 

Our unfiltered observations were published in Kiss & Derekas 
(2000). Recently, XX Cygni has been revisited by Zhou et al. 
(2002) and Blake et al. (2003). The history and observational 
record of XX Cyg can be found in those papers. Contrary to 
our earlier conclusion, these authors provided convincing ar- 
guments for slow continuous period increase with at a rate 
of about = -hi. 13 X 10"** year"' (Blake et al. 2003) and 
0.94 X lO"** year"' (Zhou et al. 2002). Here we add one night 
of observations carried out in 2002. The comparison stars were 
the same as in Kiss & Derekas (2000). 

The obtained light curves are plotted in Fig. while the 
determined epochs of maximum are shown in Table ^ The 
mean O-C value with the ephemeris from Szeidl & Mahdy 
(1981) is -t-O?0079 which falls close to any recent prediction. 
The light curve shape is similar to the ones shown in Blake et 
al. (2003) and a well-expressed post-maximum bump is present 
in the 7?-band data. Only weak evidence is present for the pos- 
sible transient secondary maximum that is discussed in Blake 
et al. (2003): there is a sHght hump around HJD 2452476.54 
which resembles that of in Fig. 2 of Blake et al. (2003) but no 
conclusive remark can be given based on this. 
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Table 5. Journal of observations and new times of maximum 
for SZ Lyn 



Date 


filter 


Inst. 


Points 




2001-02-25 


V 


P60 


395 


2451966.5244 


2001-02-26 


V 


P60 


239 


2451967.3683 


2001-04-03 


V 


Sz40 


287 


2452003.4072 


2001-05-02 


V 


Sz40 


619 


2452032.3351 



The four new times of maximum supplemented with the 
Hipparcos epoch (HJD 2448500.0560, ESA 1997) were added 
to the Hst of maxima (Paparo et al. 1988, Moffett et al. 1988). 
The resulting O-C diagram is plotted in Fig. |9l It was calcu- 
lated with the ephemeris 

HJD^^ = 2438124.39955 + 0.120534910 x E 



-1.6 



> 

<-1.4 



-1.2 



1 


1 

SZLyn 2001.05.02/03. 


/ N 
/■ 

/ 

* 


f 

\ I 

\ /■ 


1 


1 



2032.3 2032.35 2032.4 2032.45 

HJD (2450000+) 



Fig, 8. Sample Ught curve for SZ Lyn. 
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Fig. 9. The O-C diagram of SZ Lyn. The thick solid line shows 
the paraboUc part of the fit. 



3.4. SZ Lyncis 

SZ Lyncis is the best documented example of a binary pulsating 
star with very clear light-time effects. The period variations of 
the star were studied by Paparo et al. (1988) and Moffett et 
al. (1988). Since then it has been totally neglected and no new 
times of maximum have appeared in the literature, except the 
one by the Hipparcos satellite. 

We observed SZ Lyn with two telescopes on four nights 
in 2001. The exposure times were between 15 and 120 sec- 
onds, depending on the instrument and weather conditions. 
Two nearby stars were chosen as comparisons (comp = GSC 
2979-1329, V ~ 10^5, check = GSC 2979-1343, V ^ 10"8). 
More than 1500 individual points have been acquired (TablelSjl 
and a sample light curve is shown in Fig.|8] 



taken from Paparo et al. (1988). Our new observations confirm 
the consistent picture of the light-time effect thoroughly dis- 
cussed in Paparo et al. (1988) and Moffett et al. (1988). The 
parabolic and light-time fit resulted in essentially the same pa- 
rameters as those of the mentioned papers, which is not surpris- 
ing, as the overwhelming majority of the fitted data is the same. 
We could, however, improve the value of the orbital period, 
which turned out to be forb = 1 190 + 5 d (contrary to 1 177'?7 
and 1181'?1 by Paparo et al. 1988 and Moffett et al. 1988, re- 
spectively). On the other hand, we derived a slightly smaller pe- 
riod changing rate of -p^ - (+5 + 1.8) x 10"^ year"'. We note 
that the quadratic term is still too uncertain. We have experi- 
mented with different averaging of the original O-C points to 
decrease the effects of the unequal data distribution. Different 
binnings changed a lot the parabolic coefficient (almost by a 
factor of two), and that marks a certain limit in the data inter- 
pretation. 

3.5. BE Lyncis 

Kiss & Szatmary (1995) suspected binarity of BE Lyn, based 
on the seemingly cyclic O-C diagram. In order to confirm this 
preliminary result, we have continuously monitored this star 
since then. Previous observations are summarized in Rodriguez 
(1999), who did not find any long-term amplitude change of the 
light curve. 

All four instruments of this project have been used for pho- 
toelectric and CCD observations of BE Lyn. The full journal 
of observations and the list of epochs of maximum are given 
in Table|6l Because of the different detectors, various compari- 
son stars were chosen. For photoelectric photometry with Sz40, 
HIP 45515 (V = 9'J'30, spectral type F8), for CCD photometry 
with Sz40 and P60, GSC 3425-0544 (y ^ iri) were used. For 
photoelectric photometry with SNO90, HD 79763 (V = 5™95) 
served as the comparison and HD 80079 (V = 6™90) and 
HD 79439 (V - 4™82) as the check stars. Integration times var- 
ied between 8 and 120 seconds, depending on the instrument, 
detector and weather conditions. 

The 47 new epochs supplemented the earlier O-C diagram, 
and the final dataset consists of 106 points (references can be 
found in Kiss & Szatmary 1995). We plot the resulting O-C 
diagram in Fig. ^| It has been calculated with the following 
ephemeris: 



HJD„ 



2449018.2684 + 0.09586954 x E 



(Liu & Jiang 1994). Obviously, the LITE-solution of Kiss & 
Szatmary (1995) does not hold, as there is no cyclic feature in 
the present diagram. Slight variations are suggested but neither 
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Table 6. Journal of observations and new times of maximum for BE Lyn. The first five epochs were re-determined from the 
observations in Kiss & Szatmary (1995). For multicolour observations we only hst epochs based on V-band data (either from 
Johnson or Stromgren photometry). 



Date 






filter 


Inst. 


Points 


HJ fmax 


Date 






filter 


Inst. 


Points 


HJ fmax 


1995 


-01 


31 


BV 


Sz40 


267 


2449749.4651 
2449749.5622 
2449749.6562 


1999 


03 


08 


BV 


Sz40 


172 


2451246.2747 
2451246.3704 
2451246.4674 


1995 


-02 


05 


UBV 


Sz40 


48 


2449754.3545 


1999 


03 


08 


uvby 


SNO90 


62 


2451246.3702 


1995 


-02 


13 


UBV 


Sz40 


63 


2449762.4075 














2451246.4675 


1996 


-01 


31 


BV 


Sz40 


51 


2450114.5363 


1999 


04 


24 


uvby 


SNO90 


60 


2451293.3477 


1996 


-02 


25 


BV 


Sz40 


60 


2450139.4639 














2451293.4437 


1996 


-02 


26 


BV 


Sz40 


60 


2450140.4208 


1999 


04 


25 


uvby 


SNO90 


18 




1997 


-02 


20 


BV 


Sz40 


75 


2450500.5075 


1999 


05 


12 


uvby 


SNO90 


56 


2451311.4673 


1998 


-01 


30 


BV 


Sz40 


114 


2450844.3906 


1999 


05 


14 


uvby 


SNO90 


37 


2451313.4802 














2450844.4860 


2000 


02 


10 


BV 


Sz40 


57 


2451585.3663 


1999 


-01 


26 


uvby 


SNO90 


99 


2451205.5313 
2451205.6277 


2001 


01 


19 


uvby 


SNO90 


33 


2451929.4424 
2451929.5377 














2451205.7231 


2001 


02 


24 


B 


P60 


215 


2451965.2974 


1999 


-02 


27 


BV 


Sz40 


234 


2451237.2625 
2451237.3596 


2001 


02 


24 


V 


P60 


830 


2451965.3927 
2451965.4896 














2451237.4552 


2001 


03 


14 


V 


Sz40 


224 
















2451237.5500 


2001 


03 


16 


V 


Sz40 


438 


2451985.3329 


1999 


-03 


01 


BV 


Sz40 


149 


2451239.3712 














2451985.5235 














2451239.4684 


2001 


12 


08 


V 


Sz40 


193 


2452252.5191 


1999 


-03 


02 


uvby 


SNO90 


305 


2451240.4268 
2451240.5242 


2002 


04 


06 


V 


Sz28 


182 


2452371.4005 
2452371.4946 














2451240.6200 


2002 


05 


01 


V 


Sz40 


1012 


2452396.3266 


1999 


-03 


05 


BV 


Sz40 


128 


2451243.3020 
2451243.4008 














2452396.4225 



0.004 


1 1 , 1 1 1 1 1 
" C=HJD2449018.2684-I-0.09586954E 
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Fig. 10. The O-C diagram of BE Lyn 
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Fig. 11. A comparison of two simultaneous datasets of BE Lyn 
obtained with Sz40 (open circles) and SNO90 (filled circles). 



a Unear nor a parabolic fit describes very well the general ap- 
pearance. 

We have made several tests to improve our understanding 
of the limitations in the case of BE Lyn. First, we have checked 
the consistency of the determined epochs of maxima with si- 
multaneously obtained observations. On March 8, 1999, there 
were two cycles observed in parallel with Sz40 and SNO90. We 
show the light curve comparison in Fig.^ Spanish data were 
shifted vertically to the best match with the Hungarian observa- 
tions. The independent measurements yielded almost the same 
times of maximum, both were within (fiOOOl (see Table|6j. So 



we believe that the estimated accuracy of +0?0003 for the m- 
dividual epochs is indeed a realistic value. There is, however, 
a larger scatter in the O-C diagram which appears to be larger 
than the observational noise. Fourier analysis of the diagram 
did not help, because both the residuals of the linear fit and 
the parabolic fit (see below) showed peaks close to the pulsa- 
tional frequency, implying higher order period change (that is 
a simple consequence of the curved nature of the diagram) or 
low-frequency components in the spectrum caused by the inap- 
propriate fit. 
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Fig. 12. Subsequent Fourier spectra of three subsets (1995-1997, 1999 and 2001) of BE Lyn (three columns) after prewhitening 
with the dominant frequency and its harmonics. The highest peaks are marked in every spectrum, with which the next prewhiten- 
ing steps were done, from top to bottom. The small inserts show the window functions. Although there is some evidence for weak 
secondary periodicity, no features can be identified as being the same in all three subsets. The S/N ratios (Breger et al. 1993) of 
the marked peaks range from 3 to 10. 



Secondly, we have searched low-amplitude secondary peri- 
odicity in the light curves. Garrido & Rodriguez (1996) have al- 
ready tried to detect microvariability of BE Lyn but found noth- 
ing, based on five nights of observations obtained with a span of 
two years. Since we have 27 V-band light curves distributed al- 
most uniformly from 1995 to 2002, we performed their Fourier 
analysis with the same aim. We have divided the data into three 
subsets, in which the changes of the dominant frequency can be 
neglected. Nightly mean values have been subtracted (as did, 
e.g.. Kiss et al. 2002a) and that is why we had to exclude a few 
nights with short light curve coverage from the analysis. The 
first subset covered 1995-01-31-1997-02-20; the second 1999- 
01-26-1999-05-14; and the third 2001-02-24-2002-05-01. The 
most homogeneous and best-quality subset is the second one, 
which consists mainly of high-precision photoelectric photom- 
etry with SNO90. 

We prewhitened each subset with /o = 10.430842 d"' and 
its harmonics up to the fifth order and analysed the residuals for 
periodicity. The results can be summarized as follows. In every 
subset, we detected low-frequency (/ < 2d"') components 
which we attributed to nightly zero-point uncertainties. After 
their removal, we could not infer any periodicity that was com- 
mon to all subsets. The prewhitened spectra in Fig. 1121 show 
that, although there are some hints of multiperiodicity, no iden- 
tical features can be identified in all three subsets. Most impor- 
tantly, the second subset (middle column in Fig. ll2> does not 



reveal anything with amplitude larger then 5 mmag. Therefore, 
we conclude that: i) there is no detectable stable secondary pe- 
riodicity in the light curve of BE Lyn above the mmag level; ii) 
if the scatter of the 0-C diagram is real than it implies some 
light curve instabilities (transient phenomena) which are not 
strictly periodic but affect the light curve shape. Nevertheless, 
we admit that the whole effect is just on the limit of detectabil- 
ity and no unambiguous conclusion can be drawn at the present 
time. And, as noted by Garrido & Rodriguez (1996), the mi- 
crovariability can be studied only with the highest quality ob- 
servations (long-term homogeneous photometry with internal 
precision of a few mmag), which is, unfortunately, not the case 
for our data. 

We have also checked the Fourier amplitude and phase pa- 
rameters, used to describe the light curve shape, for seasonal 
variability. We have checked R21, Rn, tfii and (^31 (see, e.g., 
Poretti 2001) but found no changes. For instance, R21 was 
0.335, 0.328 and 0.303, while R^ varied from 0.108, 0.111, 
0.1 13. As a comparison, the uncertainty of these values is a few 
in the last digit and we do not consider the mentioned varia- 
tions as significant. Thus we exclude the presence of long-term 
light-curve shape changes. 

By setting aside the higher-order period change, we have 
determined an ephemeris coiTection with a linear fit and esti- 
mated the period changing rate with a parabolic fit. The linear 
approximation yields the following improved ephemeris: 
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Table 7. Journal of observations and new times of maximum 
for DY Peg 



0.005 r 



Date 


filter 


Inst. 


Points 




2001-01-01 


B 


P60 


101 


2451911.2338 










2451911.3062 


2001-08-03 


V 


Sz40 


166 


2452125.5619 


2001-08-04 


V 


Sz40 


152 


2452126.5833 


2001-08-12 


V 


Sz40 


253 


2452134.5332 










2452134.6062 


2001-08-15 


V 


Sz40 


339 


2452137.5228 










2452137.5955 


2001-09-05 


uvby 


SNO90 


57 


2452158.3805 


2001-09-08 


uvby 


SNO90 


69 


2452161.4431 



0.2 



iO.4 



> 
< 



DY Peg 2001.08.12/13. 



2134.52 2134.56 2134.6 

HJD (2450000+) 

Fig. 13. Sample light curve for DY Peg. 



HJD„ 



2449018.2681 + 0.095869521(4) x E 



which can be used to predict forthcoming maxima. The 
parabolic fit gives f3 = (-1.3 + 0.5) X 10"'^ that would corre- 
spond to - (-5 + 1. 9)xl0"** year"'. However, we consider 
this value only as an upper limit for the period change, as the 
present situation is not as clear as for the previous examples. 
Nevertheless, this limit is in good agreement with theoretical 
expectations of BP98. 



C=HJD 243275 1.9655+0.072926302 E 
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0.003 - 
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Fig. 14. The O-C diagram of DY Peg. The solid line denotes 
the parabolic fit. 



Our observations were carried out at three observatories 
on seven nights in 2001. The full log of observations, supple- 
mented with the new times of maximum, is presented in Table 
Throughout the observations we used the same compari- 
son star with all of the instruments (HD 218587, V = 9'."80, 
b - y = 0™376, mi = 0?'180, cj = 0™399, as given in 
the SIMBAD database). The secondary comparison was either 
GSC 1712-0542 (V ^ 11™7) or GSC 1712-1246 (V ~ 11™1). 
The integration time varied between 30 and 120 seconds. A 
sample light curve is shown in Fig. 1131 

After collecting and supplementing the full list of maxima 
(Blake et al. 2000 and references therein. Van Cauteren & Wils 
2000, Agerer et al. 2001) we show in Fig.UHthe final O-C 
diagram covering 1943-2001 with the ephemeris of Mahdy 
(1987): 



//7D,nax = 2432751.9655 + 0.072926302 x E 
The paraboUc fit has to following form: 



3.6. DYPegasi 

DY Peg is one of the best observed HADSs, with observations 
dating back more than a half century. A recent update of the 
period change was published by Blake et al. (2000), while the 
physical parameters of the star were determined by Wilson et 
al. (1998) and Pena et al. (1999). Garrido & Rodriguez (1996) 
presented some evidence for secondary periodicity, although 
based on only four nights of observations. The most detailed 
period study is that of Koen (1996), revealing the complex na- 
ture of the period variations. Interestingly, even the seemingly 
simple parabolic O-C descriptions yielded ambiguous results: 
Blake et al. (2000) estimated a quadratic term of the O-C di- 
agram almost a factor of seven less than, e.g. Mahdy (1987) 
(-1.36 X 10-12 vs. -0.231 X lO^'^). 



O-C = -0.0037(1) + 5.6(2) x lO^^E - 2.0(1) x lO'^E^ 

with an rms of OJOOl. This results in = (-2.79 ± 0.14) x 
10"^ year which is in good agreement with Mahdy (1987) 
and Pena et al. (1987), thus we do not confirm the different re- 
sult of Blake et al. (2000). We also note that the updated O-C 
diagram is obviously not fully described by the parabolic fit; for 
instance, the last ~ 25000 cycles would need a much steeper 
function than the present parabola. Therefore, we fully agree 
with Koen (1996) that viable quantitative alternatives to the tra- 
ditional O-C analysis must be used when investigating the fine 
details of the period changes. Here we did not go beyond the 
simplest analysis, because that would lead us too far from the 
primary purposes of this paper 
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Table 8. Journal of observations and new times of maximum 
for GW UMa 



Date 


filter 


Inst. 


Points 




2002-04-30 


V 


Sz40 


820 


2452395.5408 


2002-05-01 


V 


Sz40 


180 


2452396.5584 


2002-05-02 


Ic 


Sz40 


674 


2452397.3702 










2452397.5782 


2002-05-03 


V 


Sz40 


229 


2452398.3850: 


2002-05-03 


Ic 


Sz40 


258 


2452398.3877 


2002-05-22 


B 


Sz40 


181 


2452417.4858 


2002-05-31 


B 


Sz40 


204 


2452426.4259 



GWUMa 




I , I , I , I I I , L- 

-0.2 0.2 0.4 0.6 0.8 

phase 



Fig. 15. The phase diagrams of GW UMa (£o=2452395.5408, 
P = 0^2031940). 



3.7. GW Ursae Majoris 

The Ught variations of GW UMa ((V) « 9.7, Ay = (r44, 
P - O?2032) were discovered by the Hipparcos satellite, and 
despite its brightness and short period, there has been no new 
observation since the discovery. A special interest can be at- 
tributed to this star due to the period value: it is just on the 
boundary between HADSs and the shortest period RR Lyrae 
stars (Poretti 2001). In order to place the star in the light curve 
diagnostic diagrams used by Poretti (2001), we observed the 
star in B, V and Ic bands on 6 nights in 2002 (Poretti utilized 
I-band data of hundreds of stars to infer characteristic Fourier- 
parameters of HADS and RRc variables). More than 2500 in- 
dividual points were obtained (the full observing log is pre- 
sented in Table|8}. The comparison star was GSC 3011-2535 
(V « 101^9). 

Our analysis showed that the star is likely to be monoperi- 
odic, as we could not detect cycle-to-cycle changes of the light 
curve shape in any band (although this conclusion is based on a 
fairly scanty dataset). The phased light curves are shown in Fig. 
[Tsl The amplitudes in the instrumental BVIc system are 07^62, 
0'."44 and CPU, respectively (all but one night of observations 
were single-filtered, for better time resolution, but the trade-off 
is the lack of standard transformations). The Fourier amplitude 
parameter R21 =0.40 for the /-band data placed the star among 
the HADSs (see Fig. 7 in Poretti 2001). Therefore, we exclude 
the possibility of first- or second-overtone RR Lyrae pulsation. 



as suggested for other short-period pulsators of similar light 
curve shapes and periods (see, e.g. Kiss et al. 1999 for a dis- 
cussion of the nature of V2109 Cyg which, however, recently 
turned to be a multiperiodic HADS (Rodriguez, in prepara- 
tion)). Eight new times of maximum were determined (Table 
|S} and the mean O-C value with the Hipparcos ephemeris 
(£0 = 2448500.1160, P = 0^2031940) is -O?0062 (c.c. -9 
min), suggesting either inaccurate Hipparcos period, or slight 
period change during the last decade. Assuming a constant pe- 
riod, we calculated to following corrected ephemeris: 

//7Dmax = 2448500.1 160 + 0.20319367(3) x£ 

In any case, further CCD observations are needed to clarify 
the situation. The star is probably a Pop. II object, as the high 
galactic latitude (+59? 14) is associated with large radial veloc- 
ity (the SIMBAD database Hsts -81 km s"'). Unfortunately, 
the Hipparcos parallax is not useful (0.58 + 2.88 mas), thus no 
other information can be deduced from the presently available 
data. 

4. Summary 

In this paper, we presented observational results of an almost 
8-year long project, with the main goal of updating our knowl- 
edge on the period change of selected bright northern HADSs. 
Our results strengthen the generally adopted view of evolution- 
ary period changes, and the determined period decreasing and 
increasing rates are in good accordance with the theory. The 
most relevant new results are mostly negative results, in the 
sense of disproving earlier conclusions that appeared in the lit- 
erature: 

1 . For DY Herculis, we do not confirm the suspected binarity 
(Pocs & Szeidl 2000). The updated O-C diagram between 
1938 and 2001 does not show clear signs of a light-time 
effect. Instead of that, a quadratic fit is fairly satisfactory. 

2. For YZ Bootis, we do not confirm the increase of the period 
(Hamdy et al. 1986). Instead of that we find slow period 
decrease, close to the limits of detection. 

3. For BE Lyncis, we do not confirm the suspected bina- 
rity (Kiss & Szatmary 1995). The updated O-C diagram 
between 1986 and 2002 suggests complex period varia- 
tions with no real cyclic nature. Fourier-analysis of the 
light curves did not reveal possible multiple periodicity 
above the millimag level. A formal quadratic fit to the O-C 
diagram resulted in an upper limit to the period change 
rate that is in agreement with the theoretical expectations. 
Further observations of the star are needed. 

4. For DY Peg, we do not confirm the results on its smooth 
period change rate (Blake et al. 2000). Although there are 
hints for more complex period variation (see also Koen 
1996), the updated O-C diagram between 1943 and 2001 
can be very well fitted with a parabola that has similar co- 
efficients than in Mahdy (1987) or Pena et al. (1987). 

5. We presented the first observations of GW UMa since its 
discovery by the Hipparcos satellite. Empirical evidence 
was found for its HADS nature and the available kinematic 



A. Derekas et al.: New observations and period updates for 7 HADS 



11 



data suggest the star to be a Pop. II object (thus possibly 
belonging to SX Phe stars). The period of O'?203 makes the 
classification an intriguing task, as GW UMa is right on the 
border between RR Lyrae and HADS stars. The stabiUty of 
the light curve shape excludes the possibility of a relatively 
high-amplitude secondary period. 

For the remaining two stars (XX Cyg and SZ Lyn), our con- 
clusions are more or less confirmations of previously published 
results. As expected from the longer time basis available for the 
target stars, the newly determined quadratic coefficients (and 
the corresponding period changes) have better defined values 
and our calculations yielded a slight or significant decrease of 
the absolute value of the period changes. With these correc- 
tions, we find better agreement with theory. 

Finally, let us point out again that HADS stars are easy tar- 
gets for small astronomical instruments. The open questions 
raised for our programme stars or those of for other bright 
HADSs make continuous photometric monitoring a useful 
project for small telescopes under Ught-poUuted urban environ- 
ment. Even well-equipped amateur astronomers or small col- 
leges can make significant contribution to observational studies 
of high-amplitude 6 Scuti stars. 
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